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Abstract 



Partonic cross sections for the production of massive objects in hadronic collisions receive 
large corrections when the invariant mass of the initial-state partons is just above the 
production threshold. Since typically the center-of-mass energy of the hadronic collision 
is much higher than the mass of the heavy objects, it is not obvious that these con- 
tributions translate into large corrections to the hadronic cross section. Using a recent 
approach to threshold resummation based on effective field theory, we quantify to which 
extent the fall-off of the parton densities at high x leads to a dynamical enhancement of 
the partonic threshold region. With the example of Drell-Yan production, we study the 
emergence of an effective physical scale characterizing the soft emissions in the process. 
We derive compact analytical expressions for the resummed Drell-Yan cross section and 
rapidity distribution directly in momentum space. They are free of Landau-pole singu- 
larities and are trivially matched onto fixed-order perturbative calculations. Evaluating 
the resummed cross sections at NNNLL order and matching onto NNLO fixed-order 
calculations, we perform a detailed numerical analysis of the cross section and rapidity 
distribution in pp collisions. 



1 Introduction 



The Drell-Yan process [1], the production of a lepton pair in hadron-hadron collisions, has 
played an important role in establishing the parton picture underlying the description of hard 
interactions in QCD. In current experiments, studies of the Drell-Yan cross section as a func- 
tion of the invariant mass of the lepton pair are used to search for new heavy particles such 
as a hypothetical Z' boson, while the differential distributions provide detailed information 
about the parton distribution functions (PDFs), including in particular the sea-quark distri- 
butions. Drell-Yan production also serves as a prototype for other collider processes, such as 
Higgs production or the production of new particles. 

A lot of effort has been put in obtaining accurate theoretical predictions for the Drell-Yan 
process in perturbative QCD. The calculation of the cross section and rapidity distribution 
at next-to-leading order (NLO) in a s was accomplished in the pioneering work [2]. The first 
next-to-next-to-leading order (NNLO) result for the cross section was presented in [3] and 
confirmed much later in [4]. The rapidity distribution at NNLO was derived in [5, 6], while 
the fully differential cross section was obtained recently in [7, 8]. 

As the invariant mass M of the lepton pair approaches the center-of-mass energy of the 
collision, there is less and less phase space available for the emission of QCD radiation. After 
the cancellation of virtual and real soft divergences large Sudakov logarithms remain, because 
the scale associated with the soft radiation is much smaller than M. These "threshold log- 
arithms" threaten the convergence of the perturbative expansion and need to be resummed 
to all orders. For the inclusive cross section da/dM 2 this was accomplished in the seminal 
papers [9, 10] based on the solution of certain evolution equations in Mellin moment space 
(see also [11, 12]). The generalization of this method to the Drell-Yan rapidity distribution 
was obtained in [13]. Recent analyses of threshold resummation for the rapidity distribution 
in Drell-Yan or electroweak gauge boson production can be found in [14, 15] in next-to-leading 
logarithmic (NLL) approximation, while in [16, 17] the resummation is extended to the next- 
to-next-to-next-to-leading logarithmic (N 3 LL) order. 

Because the PDFs are strongly suppressed in the endpoint region x — > 1, the cross section 
da/dM 2 is a steeply falling function as M approaches the kinematical endpoint y/s. In fact, in 
a typical experiment it will not be possible to observe Drell-Yan pairs with masses exceeding 
about one half of the center-of-mass energy. For instance, at the LHC one does not expect 
to discover new heavy particles with masses in the 10 TeV range. In practice, one is therefore 
never in a region where the ratio r = M 2 /s approaches 1. Since threshold resummation 
deals with logarithms of the form ln(l — r), it is then not obvious why such terms should be 
treated on different footing than other higher-order terms. In view of this, it is surprising that 
large resummation effects were recently claimed to be important for the Drell-Yan rapidity 
distribution as measured by the E866/NuSea collaboration [18]. Specifically, Ref. [15] claims 
that for ^/s = 38.76 GeV and M = 8GeV (corresponding to r ps 0.04) the resummation of 
threshold logarithms at NNLL order would lower the cross section by about 30% compared with 
the fixed-order NLO result, whereas the fixed-order NNLO corrections increase it by a small 
amount [5, 6]. However, significantly smaller resummation effects were found by other authors 
for LHC kinematics (M = 115 GeV and yHi = 14 TeV) [16] and for VF-boson production at 
RHIC (M = 80.4 GeV and y/s = 500 GeV) [14]. An argument why threshold resummation 
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effects could be important even if r <C 1 has been given in Refs. [19, 20]. The idea is that 
the sharp fall-off of the parton luminosity at large x dynamically enhances the contribution 
of the partonic threshold region z = M 2 /s — > 1, i.e., the region where the center-of-mass 
energy a/§ of the initial-state partons is just sufficiently large to produce the Drell-Yan pair. 
It could then be important to resum logarithms of the form ln(l — z) in the hard partonic 
cross section. However, since (1 — z) is not related to a small ratio of external physical scales, 
it is not obvious how to give a formal justification of this argument. 

To study this question quantitatively and to assess the importance of resummation effects 
for the Drell-Yan rapidity distribution were the main motivations for the present work. To 
do so, we use a recent approach to Sudakov resummation based on effective field theory [21, 
22]. Contrary to the standard treatment in Mellin moment space, this framework completely 
separates the effects associated with different scales in the problem, thereby avoiding the 
Landau-pole ambiguities inherent in the standard approach. It then uses renormalization- 
group (RG) evolution to resum logarithms of scale ratios. The resummation is performed 
directly in momentum space, which makes it simpler to compare to and match onto fixed- 
order calculations. Our framework is particularly well suited to study the resummed rapidity 
distribution, for which we derive an exact analytic expression as a one-dimensional integral over 
PDFs. Using the convergence properties of the perturbative expansion after scale separation 
as the primary criterion, we study in detail how and under which circumstances an effective 
physical scale fi s <^ M emerges, which is associated with the soft emission in the process. Our 
approach resums logarithms of the ratio Mj \x s to all orders in perturbation theory. 

Our main findings can be summarized as follows: 

i. In the true endpoint region r — > 1, the effective soft scale \i s is an order of magnitude 
smaller than the naive guess M(l — r). For PDFs behaving like /i/jv(^) ~ (1 — x) hi 
near x — > 1, we find /j, s ~ A _1 M(1 — r) with A ~ 2 + b q + bg = O(10). This result pro- 
vides a formal justification to the argument of a dynamical enhancement of the partonic 
threshold region due to the fall-off of parton densities. 

ii. The dynamical enhancement of the threshold contributions remains effective down to 
moderate values r ~ 0.2, while at very small r values the parameter A decreases to 
about 2. This reflects the fact that for small x values the fall-off of the PDFs is much 
weaker than for large x. 

iii. Even far away from the true threshold the Drell-Yan cross section receives its dominant 
contributions from those terms in the hard partonic cross section that are leading in 
the limit z — > 1. Assuming this is true for other processes as well, the evaluation 
of virtual corrections plus soft emissions provides a simple and efficient way to obtain 
useful approximations for higher-order perturbative corrections. 

iv. With the appropriate choice of the effective soft scale fj, s , the convergence of the pertur- 
bative expansion is greatly improved by the resummation. However, for small Drell-Yan 
masses the terms beyond 0(a 2 s ) in the resummed expression for the cross section are 
numerically unimportant. We thus do not confirm the large impact of threshold resum- 
mation on the Drell-Yan rapidity distribution reported in [15]. For larger masses the 
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effects can be significant. For instance, the experiment E866/NuSea has reported data 
up to M = 16.85 GeV (corresponding to r « 0.19) [18]. We find that at M = 16 GeV 
resummation effects enhance the fixed-order predictions for the cross section by about 
25% at NLO, and 7% at NNLO. 

v. For the case of the integrated cross-section da/dM 2 , we perform a detailed comparison 
with the traditional resummation approach in moment space. Similar to the case of 
deep- inelastic scattering (DIS) in the region x — > 1, we find that the two approaches are 
equivalent up to power corrections, which turn out to be numerically small. An impor- 
tant conceptual difference is that in the effective-theory approach the running coupling 
is evaluated at physical short- distance scales depending only on the external variables s 
and M (and perhaps the rapidity Y). In this way, the Landau-pole ambiguities inherent 
in the standard approach are avoided. 

We begin our analysis discussing the structure of the hard-scattering kernels relevant for 
the Drell-Yan rapidity distribution in fixed-order perturbative QCD. In Section 3 we use 
the framework of soft-collinear effective theory (SCET) [23, 24, 25] to derive the standard 
factorization formula for the partonic cross section in the limit z — > 1 in terms of hard and soft 
functions, which we define in terms of Wilson coefficients of operators in the effective theory. 
The solutions to the RG equations obeyed by these coefficients are derived in Section 4. With 
these results at hand, we present exact analytic expressions for the resummed Drell-Yan cross 
section and rapidity distribution. A detailed numerical analysis of our results is presented 
in Section 5. After choosing the hard and soft matching scales in the effective theory by 
analyzing the perturbative expansions of the Wilson coefficient functions, we investigate the 
stability of the results under scale variations and discuss the impact of the resummation. 
Before concluding, we discuss the connection with the conventional moment-space approach. 



2 Fixed-order calculation and the threshold region 

We consider the production of a lepton pair with invariant mass M in hadron-hadron collisions 
at center-of-mass energy y/s (Drell-Yan process), focusing for simplicity on the reaction Ni + 
N 2 - ► 7* + A followed by 7*(?) — > l~ + / + . Our goal is to calculate the double differential 
cross section in the variables M 2 = q 2 and Y = \ In q ^ 9 3 , where Y denotes the rapidity of 
the lepton pair in the center-of-mass frame. Up to power corrections this cross section can 
be calculated in perturbative QCD and expressed in terms of convolutions of short- distance 
partonic cross sections with PDFs: 

dMHY = ZNM 2 s ^ J dx 1 dx 2 C ij (x u x 2: s,M : n f ) f i/Nl (x u n f ) f j/N2 (x 2 ,n f ) . (1) 

Here fi/N(x,/if) is the probability of finding a parton i with longitudinal momentum fraction 
x inside the hadron N, and ///is the factorization scale. The hard-scattering kernels CV,- 
have an expansion in powers of the strong coupling a s . The sum extends over all possible 
partonic channels contributing at a given order in this expansion. At leading order (~ a®) 
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only the channels (ij) = (qq), (qq) contribute, while at NLO (~ a s ) one must include (ij) = 

(qq), (qq), (qg), (gq), (qg), (gq) in the sum. 

It will be useful for our purposes to introduce the ratios 



M 2 

s 



M 2 



XiX 2 



(2) 



where s = X\X 2 s is the center-of-mass energy squared of the partonic subprocess that creates 
the lepton pair. This determines the maximum energy transferred to the leptons and the 
maximum invariant mass M they can have. In [5] the coefficient functions Cy are expressed 
in terms of the variable z and a second quantity 
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^e~ 2Y -z 

%2 



(1 -*)(! + ^ e -*T 



(3) 



These variables take values on the intervals < y < 1 and r < z < 1 subject to the condition 
that the parton momentum fractions 
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(4) 



do not exceed I. The allowed range for the rapidity is such that 2\Y\ < ln(l/r). We then 
define new kernels via 



C ij (x 1 ,x 2 ,s,M, [i f ) 



dz dy 



dx\ dx 2 



Qj(z,y,M, /if) 



[l-y(l-z)][l-(l-y)(l-z)] 



(5) 



At NLO the explicit results for these functions can be written in the form (with a s = a s (/j,f) 
and e g denoting the electric charges of the quarks in units of e) [5] 
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The /^-dependent terms can be derived from the fact that the cross section in (1) is scale 
independent, while the PDFs obey the DGLAP evolution equations [26, 27, 28]. The remaining 
functions follow from the symmetry relations 



Cqq — C, 



qq 5 



Cqg — C, 



qg i 



C gq — C, 



gq 



a 



qg\v->i-y 



(7) 



The hard-scattering kernels at NNLO have been calculated in [5, 6] and are available in the 
form of a computer program [29]. 

The explicit expressions for the coefficient functions given above contain terms that are 
singular in the "partonic threshold region" z — > 1, in which the center-of-mass energy of the 
parton subprocess is just large enough to create a lepton pair with invariant mass M. Indeed, 
the arguments of the logarithms in (6) suggest the relevance of two mass scales: a "hard" scale 
Hh ~ M, and a "soft" scale /j, s ~ M(l — z)/y/z = \/§ (1 — z). Physically, the hard scale is set 
by the invariant mass of the lepton pair, while the soft scale is of the order of the energy of 
the remnant jet X produced in the collision. In the region of parton kinematics where z — > 1 
these scales are separated, fih ^ Us, in which case the coefficient functions contain large 
logarithms irrespective of the choice of the factorization scale fif. Threshold resummation for 
the Drell-Yan cross section [9, 10, 11, 12] aims at resumming these logarithms to all orders in 
perturbation theory. 

Let us return to the structure of the relations (6) and identify the leading singular terms in 
the partonic threshold region. They are contained in C qq - and up to NLO multiply 5-functions 
in the variable y. Beyond NLO some of the leading singular terms in the expressions obtained 
in [5] multiply nontrivial functions of y, but since the ^/-dependence of the parton variables x\ 
and X2 in (4) is subleading in the z — > 1 limit one can always rearrange the expressions in such 
a way that the leading singular terms multiply (^-functions in y. Explicitly, we then obtain 



a 



5(y) + 5(l-y) 2rv 



qq 



e 2 q C(z,M :P i f ) + C l 



qq 



(8) 



where 



C(z,M,n f ) = 5(1 -z) + 
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+ 2 



L, 



1-z 



C F P F {z) + C A P A {z) + T F n f P f {z) 



(9) 



and we have defined L = ln(M 2 / fij) and L z = ln[M 2 (l — z) 2 / fi'jz]. The terms in the first 
line can be readily read off from (6). The two- loop coefficients Pi(z) are given in Appendix A. 
Note that the factor z in the argument of the logarithm L z could be set to 1 at leading order, 
but it is correctly reproduced by our resummation formula below and so we will keep it. The 
goal of this paper is to derive a formalism that resums these leading terms to all orders in 
perturbation theory. 
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Upon performing the integration over y, the leading singular terms in (9) give rise to the 
following contribution to the cross section: 



^thresh 47ra 2 



dM 2 dY 3N c M 2 s 

q 



^-C(z,M,ft f ) (10) 



x 



fg/Ni (Vre Y , fi f ) fg /N2 (y/7/ z e Y , fi f ) + f q/Nl {^jRj z e Y , fi f ) f q/N2 (y^e Y , fi f ) + ^ _^ 



The lower limit of the z integral is v /re =Fy , as appropriate for the two terms. At tree level, 
we recover the parton-model result 

J2 e q [/#i(^eV/)/?7A- 2 (^e" y ,/ij) + (q <-> q)] . (11) 



dM 2 dY 3N c M 2 s 

q 



We also note that integrating over rapidity, we obtain for the leading singular terms in the 
single-differential cross section 

^thresh 4?ra 2 fd Xl dx 2 



dM 2 3N c M 2 s 

q 



^24 f V 1 C ( Z ' M ' ^f) [fq/Nl( X l> N) fq/N*( X * Vf) + (<? <?)] , 

J X\ X2 



(12) 

where z = t/{x\X2), and the integration is restricted to the region where £1X2 > r. This 
result can be rewritten in the more convenient form 

^thresh 4na 2 rl fa 



dM 2 3N c M 2 s 

where 

1,1 dx 



I dz 

/ -C(z,M,Hf)f(r/z,Hf), (13) 



v—> f dx 

iVf) = 22 e2 q ~ [fq/NAx, fif) fq/N 2 (y/x, fif) + (q^ q)] (14) 

q Jy x 

denotes the Mellin convolution of the PDFs. 

In the following section we will derive a factorization formula for the coefficient C(z, M 2 , fif) 
using methods of effective field theory. The result is 

C(z, M,fi f )= H(M, fif) S(V§ (1 - z),fi f ) , (15) 

where H and S will be referred to as hard and soft functions, respectively, and will be defined 
in terms of Wilson coefficients of operators in SCET. The calculation of the components H 
and S at any order in perturbation theory is much simpler than the calculation of the Drell- 
Yan cross section at the same order. Eq. (15) thus provides an approximation to the cross 
section that requires a minimal amount of calculational work. The all-order resummation of 
the partonic threshold logarithms is then achieved by solving RG equations. 

It must be emphasized at this point that the variable z is not set by external kinematics, 
but instead is integrated over the interval between r = M 2 / s and 1. It is therefore necessary 
to specify under which conditions the partonic threshold region requires special attention. For 
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Figure 1: Fall-off of the parton luminosity function f(y, [if) for fif = 8GeV. The dashed lines 
show the asymptotic behavior for small and large y. 



where scales much smaller than the center-of-mass energy are important, such as the emission 
of a soft gluon into the final state. It is a well-known fact, however, that upon integration over 
the entire phase space the perturbation series is insensitive to scales lower than ^/s at leading 
power. 

There are two limits in which the threshold contributions are parametrically enhanced. 
First, if the invariant mass of the lepton pair is near the kinematic limit set by the total 
center-of-mass energy of the hadron-hadron collision, then r « 1 and hence z > r is always 
near 1. Threshold resummation is necessary in this case and proceeds in close analogy to the 
resummation of threshold logarithms for the DIS structure functions for x — *■ 1 [9, 10, 11, 12]. 
In practice, the region r ~ 1 is irrelevant for phenomenology, since the strong suppression of 
the PDFs near the endpoint implies a very low parton luminosity in this case. 

A second way in which the threshold contributions can be enhanced arises dynamically, 
if the weight function multiplying the hard-scattering kernel under the z-integral is steeply 
falling with (1 — z) [19, 20]. In this case threshold resummation can be justified even if r is 
much less than 1. In practice, such a behavior has to result from the fall-off of the parton 
densities with increasing x. Consider for simplicity the total cross section da/dM 2 , for which 
the relevant combination of PDFs is given by the function ff(y,[if) in (14). Figure 1 shows 
that this function is indeed very steeply falling with y. Taking [if — 8GeV, one finds that 
ff(y,fif) oc y a for y — > and ff(y,[if) oc (1 — y) b for y — > 1, where a »s —1.8 and b ~ 11. The 
figure shows that the first form reasonably well describes the behavior for y < 0.05, while 
the second form holds for y > 0.3. Using these asymptotic forms for the parton luminosity 
function, we find that for r < 0.05 



instance, also in inclusive processes such as e e 



hadrons there are kinematical situations 
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Figure 2: Comparison of the complete fixed-order results (solid lines) and the contributions 
arising from the leading singular terms (dashed lines) to the Drell-Yan rapidity distribution 
at different orders in perturbation theory. On the left we use the PDF sets MRST01LO, 
MRST04NLO, and MRST04NNLO, as appropriate to the order of the calculation; in the 
right plot MRST04NNLO is used throughout. 

while for r > 0.3 



In the first case the cross section is given by a low-order moment of the hard-scattering 
kernel C, in which case the partonic threshold region z — > 1 is not parametrically enhanced. 
In the second case, on the other hand, in the limit where we treat the exponent b as a 
large parameter, the z integral receives important contributions only from the region where 
(1 — z) < (1 — r)/b. Even for r values not near 1 there is thus a parametric enhancement 
of the partonic threshold region, which turns the threshold logarithms into logarithms of the 
exponent b. The intermediate range 0.05 < r < 0.3 is a transition region, in which the 
dynamical enhancement of the threshold region ceases to be effective as r is lowered. 

While, as we have just discussed, the dominance of the partonic threshold region cannot 
be justified parametrically for r < 0.05, it nevertheless appears that even in this case the 
Drell-Yan cross section and rapidity distribution are well approximated by keeping only the 
leading singular terms in the hard-scattering kernel (8). The reason is an inherent property of 
the hard-scattering kernel, which appears to receive the largest radiative corrections from the 
region of phase space corresponding to Born kinematics. In other words, the effects of hard 
real emissions appear to be suppressed compared with virtual corrections and soft emissions. 
To illustrate this point, we show in Figure 2 the full fixed-order predictions for the Drell-Yan 
rapidity distribution at M — 8 GeV and y/s = 38.76 GeV and compare them to the results 
obtained by keeping only the threshold terms in the hard-scattering kernel. We use the MRST 
sets of PDFs compiled in [30]. Throughout our work we use the three- loop running coupling 
a s (/x) normalized to a s (M z ) = 0.1167, and we take rif = 5 for the number of light quark 
flavors. In the right plot we use the same NNLO parton densities for all curves, so that 
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the size of the different perturbative corrections to the hard-scattering kernel can be read off 
directly. Since our main focus in this work is on the behavior of the perturbative expansion 
of the hard-scattering kernel, we will from now on always use the set MRST04NNLO for the 
PDFs. Even though r ~ 0.04 is very small in this example, we observe that at central rapidity 
about 80% of the NLO correction and 63% of the NNLO correction arise from the leading 
singular terms. For the total cross section da/dM 2 about 93% of the NLO correction and 
97% of the NNLO correction are accounted for by the leading singular terms. If we lower 
the factorization scale to Hf — M/2, then the leading singular terms come even closer to 
reproducing the fixed-order results. 

3 Derivation of the factorization formula 

The factorization theorem (15) for the leading singular terms in the hard-scattering kernel C 
has been established a long time ago [9, 10, 11, 12]. In particular, it has been understood that 
the soft function S can be represented as the vacuum expectation value of a certain Wilson 
loop of soft gluon fields. Nevertheless we find it useful to rederive this formula using methods 
of effective field theory. The advantage of this approach is that we will relate the hard and 
soft functions, H and S, to Wilson coefficients of operators in the effective theory, which obey 
certain RG equations. Solving these equations we accomplish threshold resummation directly 
in momentum space. 

In the framework of SCET the standard factorization formula (1) for the Drell-Yan cross 
section has been discussed in [31]. However, the issue of Glauber gluons [32, 33] has not yet 
been addressed in the effective theory. Glauber-gluon interactions have presented an important 
complication in the factorization analysis of the Drell-Yan process [34, 35], and it would be 
worthwhile to rederive their cancellation using SCET. As discussed below, Glauber gluons do 
not affect the factorization of the hard-scattering kernel C(z, M,/j,f) in the threshold region. 

The threshold resummation for the total cross section da/dM 2 has been performed in 
[36] using an effective-theory variant of the traditional Mellin moment-space approach. In 
the present paper we go beyond these works by presenting a derivation of the factorization 
formula (15) in SCET and performing the threshold resummation directly in z space. This 
will provide us with the tools to discuss the relevance of resummation away from the true 
threshold, i.e., for M 2 much smaller than s. 

We begin with the standard formula for the Drell-Yan cross section, 

da= ^$W)^S ( Ar i^)^ 2 (p 2 )|(-^)^ t W^(0)|iV 1 (p 1 )iV 2 (p 2 )) , (18) 

where J M = 'Y^ q e q q r y^q is the electromagnetic current. To derive the factorization theorem, 
we match the product of currents onto operators in SCET. The matching proceeds in two 
steps, and the corresponding Wilson coefficients are the hard function H and the soft function 
S, respectively. The remaining effective-theory matrix element can be identified with the 
PDFs. A similar two-step matching procedure has been used in many SCET applications. In 
particular, the factorization theorem for DIS in the region x — > 1 was derived in SCET by 
proceeding this way [22]. 
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Figure 3: Soft gluon emissions from the initial-state partons in Drell-Yan production. 



In the present paper we are solely interested in the factorization of the hard-scattering 
coefficient C(z,M,fif). For this purpose, we can simplify the effective-theory analysis by 
considering (18) with partonic instead of hadronic matrix elements. These partonic matrix 
elements correspond directly to the hard-scattering coefficients, and the hadronic cross sections 
are then obtained after convoluting the results with the PDFs. An effective-theory analysis 
of the hadronic matrix elements would be more complicated because of their sensitivity to 
nonperturbative physics governed by the scale Aqcd, which does not enter the hard-scattering 
kernels. 

Note that the product of current operators in the matrix element in (18) is not time- 
ordered. This is in contrast to the case of DIS, where the decay rate can be obtained from 
the discontinuity of a time-ordered product of currents. It is not possible to rewrite (18) in 
the same way, because the imaginary part of the time-ordered product also gets contributions 
from virtual corrections to N 1 + N 2 — > X, where X is a purely hadronic final state. A path 
integral framework to analyze operator products that are not time ordered is the Keldysh 
formalism [37, 38], which we review in Appendix C. In our discussion below, the operator 
ordering does not lead to any complications, but it implies that the soft function is not given 
by Feynman diagrams, but appropriate cuts of such diagrams, as explained in [39]. 

3.1 Derivation in SCET 

Before entering the technicalities of the discussion of factorization in SCET, it is useful to 
understand the physics behind the factorization theorem (15) by means of a simple, intuitive 
argument. Consider the special kinematics of Drell-Yan production in the (partonic) threshold 
region, as illustrated in Figure 3. Because the partonic center-of-mass energy is just above 
the invariant mass of the produced Drell-Yan pair, only soft emissions from the initial-state 
partons are allowed. As is well known, these soft emissions are described by eikonal interactions 
and exponentiate into Wilson lines. Furthermore, to leading power they leave the incoming 
partons on the mass shell, so that the production of the Drell-Yan pair is described by an 
on-shell quark form factor. 

We will now formalize this argument using the language of effective field theory. The 
effective theory is constructed by introducing fields for the momentum regions that contribute 
to the matrix elements in the given kinematics. As stressed above, for our purposes it is 
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sufficient to analyze partonic matrix elements. To distinguish the different momentum modes, 
we introduce the light-cone decomposition 

V* = (n ■ p) — + (n ■ p) — + = p\ + P- + P*i , (19) 

where n M = (1,0,0,1) and n M = (1,0,0,-1) are two light-like reference vectors (n ■ n = 
2) in the directions of the colliding partons. We denote the small expansion parameter by 
e = (1 - <C 1 and quote the components (p + ,p_,p±) of parton momenta. The relevant 
momentum regions are 

hard: p h ~ VI (1, 1, 1) , 
hard-collinear: p hc ~ VI (e, 1, y/e) , 
anti-hard-collinear: pj^ ~ VI (1, e, y/e) > 
soft: t>h * 



where the large momentum scale in the process is set by the partonic center-of-mass energy V§- 
In the first matching step, the contributions from the hard region are absorbed into Wilson 
coefficients, while the remaining contributions are represented by fields in the effective theory. 
The version of the effective theory with fields scaling in the above way is often called SCET^ 
Note that, while in many applications the soft scale is associated with nonperturbative physics, 
in our case p 2 s ~ s(l — z) 2 is of order the mass of the hadronic final-state jet and assumed to 
be in the perturbative domain. As explained above, the kinematics of the Drell-Yan process 
near threshold is such that the hadronic final state is made up of soft partons; the virtuality of 
the hard-collinear partons is parametrically larger than the jet mass: p\ c ~ p^- ~ s(l — z) ^> 
M\ ~ s(l - zf. 

The potential contribution from Glauber gluons [32, 33] presents a challenge when deriv- 
ing the standard factorization theorem of the Drell-Yan cross section into PDFs convoluted 
with hard partonic cross sections [34, 35]. However, the Glauber region does not present an 
additional difficulty for the factorization of the hard-scattering coefficient C(z, M, /if) in the 
threshold region, which we will derive below. The hard- scattering coefficient is given by the 
on-shell qq — > + X cross section, so that for its factorization only Glauber exchanges 
involving active quarks need to be studied. The momenta of the Glauber gluons relevant for 
the threshold region scale as p ~ V§(e,e, y/e), and these gluons are thus not allowed in the 
hadronic final state, which is assumed to have M\ ~ e 2 s. Writing down the loop integrals for 
Glauber exchanges between the two inital-state quarks, one finds that they are scaleless and 
vanish in dimensional regularization. In more physical terms, this means that the Glauber 
contribution is accounted for already by the standard momentum regions introduced above. 

In order not to clutter the notation, we consider for the moment a single quark flavor with 
charge e q — 1. We thus match the current J M = ip^i/j onto an effective current operator in 
SCET containing a hard-collinear quark and an anti-hard-collinear anti-quark [22, 31, 40, 41]. 
At leading power in e, this yields 

JM (0 ) _> j ds dt C v (s, t, iM f ) (ej^wy (an) ^ «^ c ) (tn) , (20) 
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where Wh c and are the usual hard-collinear (along the n-direction) and anti-hard-collinear 
(along the n-direction) Wilson lines of SCET. Note that the matching relation for the (trans- 
verse) vector current does not include a two-gluon operator at leading power in the effective 
theory. Such an operator would arise, however, in the corresponding matching relation for a 
scalar current relevant for Higgs production. Mixed quark-gluon operators corresponding to 
the (ij) = (gq), (qg), (qg), and (gq) scattering channels are power suppressed. Accordingly, 
the corresponding hard-scattering kernels do not contain terms that are singular in the 
limit z — > 1, as is explicitly seen in (6). The Fourier transform of the position-space Wilson 
coefficient 

C v (-n -pm- p 2 , n f ) = Jdsdt C v (s, t, n f ) e -^- Pl -itn- P2 (21) 

is a function of the product of the large light-cone momentum components carried by the 
quark fields. In our case, n ■ p\ n ■ p 2 = q 2 is equal to the hard scale set by the mass of the 
Drell-Yan pair. As a result the Wilson coefficient is evaluated at the time-like momentum 
transfer carried by the current, Cy(— q 2 — ie,/if). The ie prescription is required since this 
function has a branch cut along the positive real q 2 axis. Note that the same Wilson coefficient 
Cy appears in DIS, but evaluated at space-like momentum transfer [22, 36, 41]. The coefficient 
Cy can be determined by on-shell matching; indeed, it is simply given by the on-shell massless 
form factor in QCD [21]. The Drell-Yan cross section involves the current squared, so that 
the hard function in (15) is given by 

H(M,fi f ) = \C V (-M 2 - ie,fi f )\ 2 . (22) 

The expression for the coefficient Cy up to 0(a 2 ) was derived in [22, 36] and can be found in 
Appendix B. 

At leading power in e, only the n ■ A s component of the soft gluon field couples to the 
hard-collinear fields. These eikonal interactions can be represented by Wilson lines. In the 
effective theory this is achieved by redefining the hard-collinear fields as [24, 42] 

M*) - S n (x.) , Al(x) - S n (x_) A«°\x) Sl(x.) , (23) 

which implies (WU hc )(x) - S n (x__) (WU hc Y°\x). Here 

S n (x) — P exp ^ig J ds n • A s (x + sn^j (24) 

is a soft Wilson line along the n light-cone. The same redefinition, but with n and n inter- 
changed, decouples the soft gluon field also from the anti-hard-collinear fields. As a result, 
the current operator 

^(0) - (^W^f\sn)rt (StS n )(0) (wU hc f\tn) (25) 

splits into three parts, which no longer interact with each other. In the same way, the matrix 
element for the Drell-Yan process factorizes in the form 

(N 1 ( Pl )N 2 (p 2 )\(-g^)J^(x)r(0)\N 1 (p 1 )N 2 (p 2 )) 
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x(iV 2 |(^W^) (0) (0) \ (W^f\x-)\N 2 ) . (26) 

To obtain this expression, we have Fierz rearranged the fermion fields and have averaged 
over the color of the external states. We have simplified the Dirac algebra making use of the 
projection properties i/iChc — and ^Chc = Che of the hard-collinear fermion fields (and likewise 
= and i^Chc = ^hc ^ or the anti-hard-collinear fields). Also, we have neglected the power- 
suppressed dependence of the hard-collinear matrix element on x_ and x± (x + and x± for the 
anti-hard-collinear matrix element), using the fact that up to power corrections the incoming 
partons fly along the beam axis. In more technical terms, we have multi-pole expanded the 
corresponding fields to leading power [25, 43]. The soft matrix element Wuy(x, ///) (not to 
be confused with the hard-collinear Wilson lines) is a closed Wilson loop formed from the 
product of the soft Wilson lines in the two currents, 

W m {x,p f ) = ^(0|Trf [Sl(x)S R (x)] T[5j(0)5 B (0)] |0> , (27) 

where the trace is over color indices, and the operator T makes explicit that the Wilson 
lines in the complex conjugate current are anti-time-ordered. A detailed discussion of 
the operator ordering is given in Appendix C. It will become evident in Section 3.2 that the 
function W DY (x, /if) is closely related to the soft function for the Drell-Yan process. 

In the second matching step we lower the renormalization scale below the soft scale and 
integrate out the soft fields. Because the soft scale is in the short-distance domain, this simply 
amounts to a perturbative calculation of the Wilson loop (27). 

To make contact with the standard treatment, where the factorization theorem is derived 
within perturbative QCD, we note that after the decoupling transformation (23) the leading- 
power Lagrangian in each of the three sectors of the effective theory is completely equivalent to 
the QCD Lagrangian. We can thus equally well evaluate the soft function with QCD Wilson 
lines instead of soft Wilson lines. The same is true for the hard-collinear matrix elements, 
where one can replace 

(^i|(a c ^ c ) ( °H^)|(<a c ) {0) (0)|iVi) - (^Vi|^ + )| [* + ,0]#))|Aq>, (28) 

where tp(x) is the QCD quark field and [x + ,0] a QCD Wilson line of n • A gluons extending 
from to x + . 

We have simplified the effective-theory treatment by restricting ourselves to on-shell par- 
tonic instead of hadronic matrix elements, which has the advantage that we do not need to 
discuss momentum regions that scale with powers of Aqcd. 1 These momentum regions cannot 
enter the hard or the soft functions, and given that the Drell-Yan cross section is known to 
factorize for s(l — z) ^> Aq CD , we are guaranteed that their presence will not spoil factoriza- 
tion. The restriction to partonic matrix elements can thus safely be dropped. The hadronic 

1 For DIS at x — > 1, a full factorization analysis (including messenger fields [42]) was performed in [22]. 
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matrix elements of the form (28) are the usual PDFs [44] 

f q/N (x, »,) = ± J°°dte- ixt ^ (N(p)\i>(tn) | [in, 0] ^(0) \N(p)) . (29) 

Inserting the factorized form of the matrix element (26) into expression (18) for the cross 
section, we obtain 

da 4-7ra 



dM 2 3N r M 2 



x 



C v (-M 2 -?e,/i/)| 2 ^ej / dx x dx 2 [f q / Nl ( x i, A*/) fq/N 2 ( x 2, ^/) + (? ^ ?)] 
/ (2^ 27 / ^ W^s) . (30) 



From now on we will drop the ie prescription in the argument of Cy. 

3.2 Kinematic simplifications in the threshold region 

To proceed, we need to study the kinematics in the partonic threshold region, where z — > 1. 
With the expansion parameter e = (1 — z) <C 1, it is easy to show that in the parton center- 
of-mass system 



151 <-y(l-*) = O(e) =► q° = ^s + 0(e). (31) 

It follows that g° is parametrically larger than the spatial components of q^ . This observation 
implies that the term q° in the denominator in (30) is independent of q at leading power, in 
which case performing the integral over d 3 q yields S^ 3 \x). It also follows that the rapidity of 
the Drell-Yan pair in the parton center-of-mass system vanishes at leading power, 

1 o° + o 3 

Y CMS = -ln\^L = O(e) , (32) 

2 q u — q 6 

which explains why the rapidity in the laboratory frame, Y = Y C ms + f m i s determined 
by the ratio x±/x 2 at leading power in (1 — z), as can be seen from (4). Finally, we need that 



(x lPl . + x 2 p 2+ - q)° = (1 - z) + 0(e 2 ) . (33) 
Using these results, we obtain at leading power in the z — > 1 limit the expression 

\C V (-M 2 , fi f ) | 2 e \ J dx ± dx ? [fq/Ni (xi, Hf) fq/N 2 (x 2 , Hf) + (q <-> q)] 



dM 2 3N r M 2 



x^ f |^e^ 1 -^°/ 2 ^ DY (a; ,f=0, / . / ), (34) 
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where the integration region over the parton momentum fractions is such that X\X 2 > r, so 
that z = rjx\x 2 < 1. In the final step we introduce the Fourier transform of the position-space 
Wilson loop at time-like separation via 

W BY (u,fi f ) = J ^e^°/ 2 W DY (x°,x = 0,fi f ). (35) 

This Wilson loop plays the role of the jet function in DIS in the sense that it describes the 
properties of the hadronic final state. It has been introduced previously in [12, 39]. The 
Drell-Yan cross section now takes the form 

da Ana 2 . , 2 ,,,v^ 2 f dxi dx 2 



dM 2 3N c M 2 s 



/dec ^ diijC^ 
— — [fq/Ni (zi, fif) fq/N 2 (x 2 , Hf) + {q^ q)] 
X\ x 2 



xVsW DY (Vs(l- z),fi f ), (36) 

and from comparison with (12) and (15) we identify the soft function as 

S(VI (l-z),n f ) = Vl W BY (V§ (l-z),Hf). (37) 

We should mention that at leading power in (1 — z) the argument of the soft function could 
be simplified as y/I (1 — z) — M(l — z)/y/z ~ M(l — z); however, since the exact expressions 
(6) for the hard-scattering kernels at NLO contain the logarithm L z = ln(s(l — z) 2 /^ 2 ), we 
prefer to keep the argument of the soft function in the form written above. 



4 Momentum-space resummation at large z 

At this point we have identified the two components in the factorized expression (15) for the 
hard-scattering coefficient C(z, M,fj,f) with field-theoretic objects defined in terms of operator 
matrix elements. The resummation of threshold logarithms arising in the z — > 1 region can 
now be accomplished by solving the RG evolution equations obeyed by these quantities. 



4.1 Evolution of the hard function 



The evolution equation for the hard matching coefficient Cy evaluated at time-like momentum 
transfer and its solution can be obtained from the corresponding results valid for space-like 
momentum transfer [21] by analytic continuation. This leads to 



d 



dlnfj, 



C v {-M 2 -ie,n) = 



cusp 



C v {-M 2 -ie,n) . (38) 



We have reinserted the ie regulator, which determines the sign of the imaginary part of the 
anomalous dimension. The appearance of the logarithm and its coefficient, the cusp anomalous 
dimension r cusp [45, 46], can be explained using arguments presented in [40]. This term in the 
evolution equation is associated with Sudakov double logarithms. The remaining term, j v , 
accounts for single-logarithmic evolution. 



15 



The exact solution to (38) is 

/ M 2 \ ~ av ^ h ^^ 

C v (-M 2 ,/j, f ) = exp [2S(fi h ,tJ>f) ~ CLyv(fih,fJif) + inarch, A*/)] ( — C v (-M 2 ,/j, h ) , 

(39) 

where \ih ~ M is a hard matching scale, at which the value of Cy is calculated using fixed- 
order perturbation theory. Note that the Wilson coefficient at time-like momentum transfer 
is a complex quantity. The Sudakov exponent S and the exponents a n are given by [47] 



ct s (p) a oc s (fi) 



S(u, /i) = - 




a a {y) a s [y) a s (v) 



and similarly for the function a 7 v. The perturbative expansions of the anomalous dimensions 
and the resulting expressions for the evolution functions valid at NNLO in RG-improved 
perturbation theory are collected in Appendix B. 



4.2 Evolution of the soft function 

In order to derive the evolution equation for the soft function it is important to have a 
consistent definition of the threshold region. To this end we consider the limit M 2 — > s 
in (36), which implies z — > 1, since z > r. The condition xix 2 > r then implies that x± t2 — > 1, 
in which case the DGLAP evolution for the PDFs can be written in the simplified form 

d , ( \ a, v , , \ . riT-i / \ [ l dz f q/N (x/z,fj) 



f q/N (x, = 2^(a s ) f q/N (x, /x) + 2r cusp (a s ) / - ^ x ' z ^ + . . . , (41) 



<iln/i 

which is obtained by expanding the Altarelli-Parisi splitting function as 

*W*) = + 27^(« s ) 5(1 -,) + ... . (42) 

I 1 - + 



This asymptotic form of the splitting function holds to all orders in perturbation theory 
The cusp anomalous dimension and the coefficient 7* have been calculated at three-loop order 
[49]. They are given in Appendix B. 

Requiring that the Drell-Yan cross section in the threshold region be independent of the 
arbitrary factorization scale /if, we find that the momentum-space Wilson loop obeys the 
integro-differential evolution equation 



dWnY(v, fx) 



AT cnsp (a s ) ln- + 2 7 w (a 



dln/i 

- ir, 



where 



W dy (uj,ij,) 

(a.) r du/ ^M)-w m (^) < (43) 

Jo to — to 



cusp 
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Curiously, this quantity starts only at two-loop order. The analogous relation in the case of 
DIS is 7 J = 7^ + 7 y , where the anomalous dimension 7 J of the DIS jet function is defined 
in analogy with ^ w in (43). Combining the two relations we obtain 2 / ~f J — 7^ = 7 y , which 
is an exact relation that links the anomalous dimensions of the soft function in Drell-Yan 
production and the DIS jet function to the anomalous dimension of the vector current. 
The exact solution to the evolution equation (43) can be written in the form [21, 50] 

1 / uo \ 2v e~ 2lE7) 

W DY (uj,iif) = exp [-AS(/i s , /if) + 2a jW (ii s , /if)] s DY (d v ,/i s ) - f-j — — , (45) 

where d v denotes a derivative with respect to an auxiliary parameter r], which is then identified 
with 77 = 2ar(/x s , ///). This result is well defined for rj > 0. The solution for negative 77 is 
obtained by analytic continuation. For instance, to obtain the result for — | < 77 < we use 
the identity 

r^iM = ^/M-/(0) + M^, (46) 
Jo oo l - 2r > Jo ^- 2r > 2?7 y 1 

where f{oj) is a smooth test function. For 77 < — | additional subtractions are required. 

The function Sdy is obtained from the momentum-space Wilson loop by the Laplace trans- 
formation [21] 

f°° 1 
s DY (L,fi s )= due S0J W DY (uj,fx s ) , s = — E7 -. (47) 

The position-space Wilson loop on the right-hand side of (35) can be shown to have the 
functional form ^^(a; ,^ = 0, /i s ) = /(|x°/i s e 7B , a s (/i s )) [12, 39], and a straightforward 
calculation shows that the function s D y(£, /-t/) can be expressed in terms of f(t,a s ) as 

s DY (L,/i s ) = /(e- L / 2 ,a s (/i s )). (48) 

We can then use the explicit two-loop expression for the position-space Wilson loop obtained 
in [39] to compute the function 3dy at two-loop order. The result is presented in Appendix B. 

4.3 Resummation of large logarithms 

We are now in a position to derive the RG-resummed expression for the hard-scattering coef- 
ficient C(z, M, /if) in (15), which is given by the product of the solutions (39) and (45) for the 
hard and soft functions. The result can be simplified by eliminating the anomalous dimension 
•y w using (44), and combining the Sudakov exponents using the relation 

S(fi h ,fif) - S(fi s ,fi f ) = S{n h ,n s ) - a T (n 8 ,iif) In — . (49) 

We find 

C(z,M,Hf) = \C v {-M 2 ,fi h )\ 2 U{M,^ h ,^ s ,fi f ) 
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Table 1: Different approximation schemes for the evaluation of the resummed 
cross-section formulae 



RG-impr. PT 


Log. approx. 


Accuracy ~ a n s L k 


r 

1 cusp 




Cy, sdy 




LL 


k = 2n 


1-loop 


tree-level 


tree-level 


LO 


NLL 


2n - 1 < k < 2n 


2-loop 


1-loop 


tree-level 


NLO 


NNLL 


2n - 3 < k < 2n 


3-loop 


2-loop 


1-loop 


NNLO 


NNNLL 


2n - 5 < k < 2n 


4-loop 


3-loop 


2-loop 



where r] = 2ar(// s , ///), and we have defined the evolution function 

— J exp [4S(fi h ,fi s ) -2(L r v(fi h ,fi s ) + 4a^(fi a , fif)] . (51) 

J 

As before, equation (50) is valid for rj > (/j, s >///)■ For negative rj (/if > // s ), integrals of 
ln n (l — z)/(l — z) l ~ 2r] with test functions f(z) must be defined using a subtraction at z — 1 
and analytic continuation in 77. 

We emphasize that the result (50) is formally independent of the scales fih and fj, s , at which 
the matching conditions for the hard and soft functions are evaluated. On the other hand, 
the hard-scattering kernel C does depend on the factorization scale at which the PDFs 
are renormalized. In practice, a residual dependence on the matching scales arises when the 
perturbative expansions of the matching coefficients and anomalous dimensions are truncated, 
and this dependence can be used to estimate the remaining perturbative uncertainties. Setting 
the three scales /i s , and /if equal to each other in the resummed expression (50), one can 
readily reproduce the leading singular terms for z — > 1 in the fixed-order perturbative QCD 
expression for the hard-scattering kernel. In this way we have obtained the two-loop corrections 
in (9). 

The final expression (50) for the hard-scattering kernel can be evaluated at any desired 
order in resummed perturbation theory. Table 1 shows what is required to obtain different 
levels of accuracy. In this work we adopt the counting scheme of RG-improved perturbation 
theory, where at LO one includes all 0(1) terms, at NLO one includes all 0(a s ) terms, etc. The 
large logarithm ln(fih/fi s ) is counted like 0(l/a s ). In the literature on threshold resummation 
the alternative notation N n+1 LL is often used instead of N"LO. The leading logarithmic (LL) 
approximation is listed only for completeness, as it misses some 0(1) terms. 

In the following section we will perform a detailed numerical analysis of the Drell-Yan cross 
section and rapidity distribution. In most cases of phenomenological relevance the invariant 
mass of the Drell-Yan pair will be small compared with the center-of-mass energy, i.e. r = 
M 2 1 s 1. Nevertheless, it is interesting to briefly consider the limit r — > 1, in which the need 
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for threshold resummation is justified parametrically (see the discussion in Section 2). In this 
case the convolution integrals in formula (12) for the Drell-Yan cross section can be performed 
analytically if a reasonably simple model for the PDFs near the endpoint is adopted. We 
parameterize the behavior near x = 1 as 



(52) 



and similarly for the anti-quark distribution. It then follows that at leading power in (1 — y) 
the parton luminosity function defined in (14) is given by 



= »M P " » 



(53) 



Furthermore, the sum over flavors is dominated by the terms with the smallest exponent 
(b q + bq). For example, in pp collisions at one finds b u + b a xs 14.4, b d + b d ~ 10.8, and 
6 S + 6s « 16.0 at /x = 3GeV [51]. These exponents increase by a flavor-independent amount 
2ctr(/i/,/io) when /iy is raised to larger values, e.g., by about 0.4 at fif = 8GeV [22]. It follows 
that the leading behavior near the endpoint is due to the down-quark contribution. 
Introducing the Drell-Yan i^-factor as the ratio 



da 



K(M 2 ,t) 



da 



dM 2 



LO 



we then obtain at leading power 



K{M\t) = \C V (-M 2 , (j, h ) | 2 U (M, n h , [j, s , fj,f) 

x (1 - r) 2 " s DY ( In M2(1 2 ~ T)2 + d v , /i. 



e-a™ r(2 + b d + b s ) 



(54) 



(55) 



T{2 + b d + b d + 2r 1 ) 

This may be compared with the i-T-factor for DIS at large Bjorken x and momentum transfer 
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Q 2 = —q 2 , which has been derived in [21] and reads 



K(Q 2 ,x) = \C V {Q 2 , n h )\ 2 U{Q, fi h , fj, a , fif) exp [-2a 7 <^(// s , ///)] 

x (1 - x)" JDIS ( In + 9„ T{l + bu + r]) , (56) 

where in this case the dominant contribution comes from the valence up-quark, which has 
the smallest b q parameter, b u ~ 4.0 at fif = 3 GeV [51]. Obviously the structure of threshold 
logarithms is very similar in the two cases once we consider equal hard scales (Q 2 = M 2 ) 
and compare the small parameter (1 — x) in DIS near the endpoint with the small parameter 
(1 — r) 2 in Drell-Yan production. Differences arise from the following facts: (i) Drell-Yan 
production has time-like kinematics, whereas DIS probes the nucleon at space-like momentum 
transfer. This gives rise to a difference in the hard matching coefficients \Cy\ 2 starting at 
one-loop order, (ii) The soft/jet functions in the two cases are different starting from one-loop 
order, (iii) The Drell-Yan cross section involves a convolution with two PDFs, whereas in DIS 
a single parton density appears. This explains the different coefficients in front of and r\ in 
the expressions for the .fT-factors. (iv) For the same reason, the resulting convolution integrals 
over the PDFs give rise to different expressions involving the b q and bq exponents. 

There is one more important piece of information that we can extract from the result (55) 
for the Drell-Yan X-factor. As we have seen, the exponents b q and bq take rather large values. 
Therefore, the arguments of the T-functions in (55) contain the large quantity (2+b d +b d ) ~ 13. 
It is straightforward to show that the derivative with respect to rj in the argument of the soft 
function has the effect of changing the argument of the logarithm as follows: 

, M 2 (l-r) 2 , M 2 (l-r) 2 

In — — + d„ -> In — — ^ J —- , (57) 

A " ^(2 + b d + b d ) 2 1 J 

up to 0(1) factors. It follows that a proper choice for the soft matching scale near r — > 1 is 

a ~ M(1 ~ r) ~ M ^~ T) (58) 
(2 + 6d + 6j) la » ( 58 ) 

which is an order of magnitude less than the naive choice M(l — r). The fact that the fall-off of 
the parton densities strongly favors the large-z region leads to a strong additional suppression 
of the effective soft scale. 



5 Phenomenological analysis 

In this section we perform a detailed numerical analysis of our results. One of our goals is 
to study to what extent threshold resummation is important (or even justified) in processes 
where the invariant mass of the Drell-Yan pair is not very close to the center-of-mass energy. 
We have seen in Section 2 that for very small values of the ratio r = M 2 /s the threshold 
contributions are not parametrically enhanced. Even though empirically these terms still 
give rise to the dominant contributions to the cross section, there is no need to perform a 



20 



3 
2 
1 

e 

u 

-1 

-2 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Hh/M 

Figure 4: Dependence of the first two expansion coefficients in the perturbative series for the 
hard function on the matching scale Hh- 

resummation of the threshold logarithms. On the other hand, the result (58) derived in the 
previous section shows that near the true endpoint the threshold logarithms are enhanced by 
two effects: the kinematic restriction of the ^-integral to the small interval between r and 1, 
as well as the dynamical enhancement from the strong fall-off of the parton densities. As a 
result, the appropriate value of the soft matching scale is an order of magnitude smaller than 
the naive choice M(l — r). From our discussion so far it is not obvious how to interpolate 
between these two extreme cases. We now perform a detailed numerical study to assess the 
importance of resummation at intermediate values of r. 

5.1 Choices of the matching scales 

We begin with a discussion of the proper choice of the matching scales Hh and fi s , using the 
convergence of the perturbative expansions of the matching coefficients Cy and 3dy i n the re- 
summed hard-scattering kernel (10) as the primary guiding principle. While it is obvious that 
the hard matching scale should be chosen of order M, the choice of the soft scale is more prob- 
lematic. Naively, based on the structure of the result (50) one would expect that ln(/i^/M 2 ) 
should in some sense be identified with the "average" value of ln[(l — z) 2 /z}. Unfortunately, 
however, the distribution in the variable z is both singular at z — 1 and not positive definite, 
so that it does not lend itself to a probabilistic interpretation. A simple way to avoid large 
logarithmic contributions would be to make the scale choice n 2 /M 2 oc [(1 — z) 2 /z] inside the 
^-integral. However, this would lead to Landau-pole singularities in the integrand and hence 
would upset the proper scale separation that is at the heart of our approach. 

Our approach to threshold resummation based on effective field-theory methods applied 
directly in momentum space provides a natural way of resolving this question. The matching 
scales Hh and fi s should be chosen such that the perturbative expansions of the Wilson coef- 
ficient functions Cy and Sdy i n (50) are well behaved. For the case of the soft function, this 
criterion should be applied after the integration over z in (10) has been performed. This is 
the essence of effective field theory: one performs matching calculations at scales where these 
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Figure 5: Relative contributions to the Drell-Yan cross section da/dM 2 at M = 20 GeV arising 
from the one- loop (left) and two-loop (right) corrections to the soft function s DY , as a function 
of the soft matching scale fi s . The curves are labeled by the corresponding values of r = M 2 /s. 



calculations can be done in fixed-order perturbation theory, and use the renormalization group 
to perform the evolution ( "running" ) between the different matching scales. 
We begin by applying this criterion to the hard function 

oo 

H(M,fi h ) = \C V (-M 2 ,fi h )\ 2 = 1 + J2c n (^) kWf. (59) 

71=1 

Figure 4 shows the dependence of the expansion coefficients c\ and C2 on the ratio fih/M. The 
one-loop coefficient c\ vanishes for Hh/M ~ 1.569 and ^h/M « 0.142. The second solution 
is in a region of very small fih, where the expansion coefficients vary strongly and where the 
ln 2 (M 2 //-4) and ln(M 2 //i^) terms have opposite sign. We will thus discard it. In the region 
around the first solution the two-loop coefficient c<i is stable and positive. In our numerical 
analysis we will vary fi^ between M and 2M, taking ^ = 3M/2 as the default choice. 

The matching scale /x s must be determined separately for each process, since it is sensitive 
to the integration range of the z variable (which depends on r and Y) and to the shape of the 
PDFs. In Figure 5, we plot the relative contributions to the cross section da/dM 2 (normalized 
to the total cross section) arising from the one- and two-loop terms in the soft function Sdy as 
a function of /j, s . We choose M = 20 GeV and consider different values of r = M 2 / s between 
0.01 and 0.7. The plots have been obtained by setting the factorization scale equal to M and 
using MRST2004NNLO parton densities [30]. We have checked that virtually indistinguishable 
results are obtained when /// is varied by a factor of 2. Notice that with increasing r values 
the regions where the one- and two-loop contributions are of modest size shift toward lower 
/i s values. To be specific, we consider two criteria for a good convergence of the perturbative 
expansion (see the left plot in the figure): 

I. Starting from a high scale, we determine the value of fi s at which the one-loop correction 
drops below 15%. 

II. We choose the value of fi s for which the one-loop contribution is minimal. 
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Figure 6: Scale-setting results for the soft matching scale fi s for different values of r = M 2 /s 
and Drell-Yan masses M = 100, 50, 20, 10 GeV. The lighter curves correspond to lower scales. 
The upper set of curves corresponds to convergence criterion I, the lower one to criterion II. 



Note that with either choice the two-loop corrections at the corresponding fi s values are very 
small, indicating that the first two terms in the perturbation series for the soft function are 
well behaved for the same choice of scale. The same analysis can be repeated for different 
masses of the Drell-Yan pair. 

The resulting values for the soft scale fi s determined using these convergence criteria are 
shown in Figure 6, where we consider the choices M = 10, 20, 50, 100 GeV. To a good 
approximation the curves for the ratio /x s /M as a function of r exhibit scaling, i.e., they are 
almost independent of M. We thus obtain a relation of the form fi s /M = g(r). Small scaling 
violations arise from the scale dependence of the PDFs and of the running coupling in the 
perturbative expansion of the soft function. Our numerical results are well reproduced by the 
empirical functions 

t Mil - t) , TT Mil - t) 

for the two criteria. Their form should not be taken too seriously except to note that for r — ► 1 
both functions approach fi s = const, x M(l — r), as required by the resummation formula (55) 
for the .ff -factor valid near the true endpoint. Indeed, the smallness of the constants in the 
two cases (0.125 and 0.080, respectively) is in good agreement with our estimate in (58). It 
results from the dynamical suppression provided by the PDFs. In the opposite limit r — ► 0, 
both forms yield fi s = const, x M with an 0(1) constant, as required by the fact that there is 
only a single physical scale in this case. Below we will vary the soft scale between f/ s and /x* 1 
and, somewhat arbitrarily, use the average of the two scales as the default choice. 

The same analysis can be carried out for the double differential decay rate cPa/dMdY , and 
it leads to similar results. In general, one finds that at central rapidity (Y = 0) the resulting 
values for the soft matching scale are higher by about 20% than in the case of da/dM 2 , while 
at larger rapidity they are lower, so that on average, after integration over Y, the results for 
the total cross section are reproduced. 
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5.2 Scale dependence and impact of resummation 

We now proceed to study the stability of our resummed expression for the Drell-Yan cross 
section, using scale dependence as an estimator of yet unknown higher-order perturbative 
effects. We vary the matching scales /ih and /i s about the default values \ih = 1.5M and 

= (/4 + / i s I )/2 determined in the previous section. To simplify comparisons with the 
literature we adopt the conventional choice /if — M for the factorization scale, at which the 
PDFs are renormalized. From the point of view of effective field theory it would be more 
natural to choose a lower value for /if, given that the cross section is sensitive to physics at 
scales much below the hard scale M. We will see, however, that our results are very stable 
with respect to variations of the factorization scale, so that the choice of the default value 
is not particularly important. Since our focus is on the behavior of different perturbative 
approximations to the hard-scattering kernel, we use the same set of PDFs (MRST04NNLO) 
throughout the analysis. 

Our results are shown in the first three plots in Figure 7, in which we study the Drell-Yan 
.fT-factor defined in (54) for M = 20 GeV and various values of r = M 2 /s. In the calculation of 
the i^-factor we keep the factorization scale Hf — M fixed in the leading-order expression in the 
denominator, even when /if is varied in the numerator. For the time being we only include the 
leading terms in the z — > 1 limit, corresponding to the result (50). Adding the small power- 
suppressed corrections would not change any of our conclusions. We observe an excellent 
convergence of the perturbative expansion for the cross section after resummation. The bands 
corresponding to the LO, NLO, and NNLO approximations overlap, and the dependence on 
the matching scales jih and ji s becomes negligible beyond LO, indicating that the residual 
perturbative uncertainty is very small. 

The third plot, showing the dependence on the factorization scale /i/, requires some com- 
ments. We first note that in the approximation where we resum the leading singular terms 
in the partonic cross section near z — > 1 but neglect power-suppressed terms, our results for 
the cross section are no longer strictly independent of the factorization scale /if. In order to 
make them formally scale invariant one should add back the power-suppressed terms in fixed- 
order perturbation theory. Despite this, we find that the resummation greatly reduces the /// 
dependence compared with the fixed-order calculation, for which results are shown in the last 
plot in the figure. (The fixed-order results can be obtained by setting /i h — /i s — /if in the 
resummed expression.) The reason is that already the LO result after resummation compen- 
sates the leading scale dependence of the PDFs through the /if dependence of the functions 
a 7 «(/i s , /if) and rj = 2ar(/J, s , fJ>f) in (50). Indeed, we see that for modest values < r < 0.2 
the dependence on the factorization scale after resummation is almost absent already at LO. 
On the other hand, there is a sizable scale dependence in the fixed-order calculation even 
beyond LO. We emphasize that there is very little experimental information about the rele- 
vant parton densities (the sea-quark distributions, in particular) at large x values, and as a 
consequence one cannot trust the numerical values of the if-factor at values r > 0.2. Indeed, 
we find that the results for the cross section at r = 0.3 obtained using different sets of PDFs 
(MRST04NNLO [30], MRST01NNLO [52], and CTEQ6.5 [53]) differ by a factor 4. 

In Figure 2, we observed that the Drell-Yan cross section is dominated by the singular 
threshold terms. We now assess whether the threshold contribution contains large logarithms 
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Figure 7: Dependence of the resummed Drell-Yan cross section for M = 20 GeV on the 
scales fih, Us, and fif. The bands show the i^-factor obtained at LO (light), NLO (medium), 
and NNLO (dark). The last plot shows for comparison the /i/ dependence in fixed-order 
perturbation theory. 



which should be resummed. In Figure 8, we compare the result for the total cross section 
obtained with the default values of the matching scales fih and /i s to the evaluation of the 
threshold terms in fixed-order perturbation theory. The differences between the two sets of 
curves show the effect of the resummation. For illustration purposes we consider two examples. 
The first is the case of pp collisions at yfs = 38.76 GeV, corresponding to the energy of the 
fixed-target experiment E866/NuSea [18]. In this experiment Drell-Yan masses in the ranges 
4.2-8.7 GeV and 10.85-16.85 GeV have been observed. As a second example we consider Drell- 
Yan production via a virtual photon at the LHC energy a/s = 14 TeV (including the Z° channel 
would not alter our results for the i^-factor significantly). The figure shows that resummation 
accelerates the convergence of the perturbative expansion. On the other hand, the plots also 
show that the most important logarithmic corrections are contained in the fixed-order NNLO 
results, at least for moderate lepton-pair masses. For large masses the resummation becomes 
more important. In the first example, the resummed cross section is about 8% (28%) larger 
than the fixed-order result at NNLO (NLO) for M = 16.85 GeV. At the LHC, on the other 
hand, resummation effects beyond the NNLO fixed-order result remain small even at very 
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Figure 8: Resummed (solid lines) versus fixed-order results (dashed lines) for the i^-factor as 
a function of M. The light, medium, and dark lines correspond to LO, NLO, and NNLO, 
respectively. Default values are used for all scales. 



large values of M. 

5.3 Rapidity distribution and cross section at y/s = 38.76 GeV 

As a final application, we now return to the rapidity distribution in Drell-Yan production at 
y/s = 38.76 GeV. As mentioned in the Introduction, in this case large resummation effects 
were found for M = 8 GeV [15] even though r « 0.04 is very small. These effects were claimed 
to reduce the NLO fixed-order cross section by about 30%. Fixed-order predictions for the 
rapidity distribution up to NNLO were discussed in [5, 6]. Here we present results for the two 
cases M = 8 and 16 GeV. In order to obtain the best possible predictions we combine our 
resummed result for the cross section with the power-suppressed terms calculated in fixed-order 
perturbation theory. In our approach this matching can be implemented in a straightforward 
way as follows: 

^0-combined ^^.thresh 

dM 2 dY = dM 2 dY 

In Figure 9, we compare our RG-improved results with the fixed-order results, varying the 
scales over the ranges M/2 < /// < 2M, M < ji h < 2M, and // s < ji s < /i* 1 . The bands 
reflect the variations about the default value. In the fixed-order case only the first variation 
is relevant, while in the resummed case we add the individual variations in quadrature. 

We observe again that resummation significantly accelerates the convergence of the per- 
turbative expansion. Moreover, even though in the resummed case we include the scale de- 
pendence from the variation of three different scales, the combined uncertainty at NLO and 
NNLO is significantly smaller than in the fixed-order case. Also, given the better overlap of 
the bands in the resummed case, our error estimates appear to be more conservative. As a 
final comment, we note that for M = 8 GeV the resummed results at NLO and NNLO are 
consistent within errors with the fixed-order results, indicating that threshold resummation 
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Figure 9: Fixed-order (Y < 0) versus resummed (Y > 0) predictions for the rapidity distribu- 
tion at y/s = 38.76 GeV and two values of M, at different orders in perturbation theory. The 
bands reflect the combined scale dependence. LO bands are light, NLO bands are medium, 
NNLO bands are dark. 



is not an important effect. This is in stark contrast to the conclusion reached in [15]. For 
the higher mass M = 16 GeV, the two NNLO bands are consistent with each other at central 
rapidity, but the resummed result is significantly higher than the fixed-order prediction for 
Y > 0.3. For the integrated cross section at this value of M, threshold resummation enhances 
the fixed-order value by about 7%. This can be seen from Table 2, which shows our final 
predictions for the integrated cross section da/dM 2 . Besides the results obtained with and 
without resummation, we also give the contributions of the resummed threshold terms alone, 
corresponding to the first term in (61). 



5.4 Resummation in moment space 

Traditionally, resummation is performed in moment rather than momentum space [9, 10]. For 
the Drell-Yan cross section integrated over rapidity one takes moments in r at fixed M: 

For the moment-space analysis of the rapidity distribution one performs a Fourier transform 
in the rapidity in addition to taking moments in r [13, 15]. In the following, we will restrict 
ourselves to the integrated cross section for simplicity. Using the representation (12), the cross 
section in moment space factorizes as 

°» = E < [/#£ f$$ + (*~ 9)] C N+1 (M\ N ) , (63) 

c q 

where the moments of the hard-scattering coefficient and the PDFs are defined in analogy 
with (62). In order to accomplish the resummation for the moments of the hard-scattering 
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Table 2: Predictions for the Drell-Yan cross section dajdM 2 for y^s = 
38.76 GeV and two values of the invariant mass M of the lepton pair. Units 
are pb/GeV 2 and fb/GeV 2 , respectively. 



M Method LO NLO NNLO 

8 GeV combined 0.436+°$ 2 0.493 0.512t°;°° 2 

threshold 0.436 t^ 2 0.482^ 0.501 +° Z 

fixed-order 0.299+°;°^ 0.449 0.505 ^11 

16 GeV combined 1.49 1.61 ±g;gj l-68±g;gj 

threshold 1.49 159^ 1.63+°;°} 

fixed-order 0.76 + °; 2 ^ 1-29^^ 1.57l°;°* 



coefficient, we go back to expression (15) and use that [54] 

^ / M 2 

dzz N - 1 S(M(l -z),n) = s DY [ ln- 

'o 



/ 1 / M 2 \ / 1 \ 

cfc /- 1 S(M(1 -z),») = s DY ( In p)+0 (-) , 



(64) 



with N = e lE N. We then insert the solution to the RG evolution for the function s D y, which 
as shown in [22] obeys a RG equation analogous to (38). This leads to 



/ M 2 \ 

Cn(M 2 , fXf) = \Cv{—M 2 ,^ h )\ 2 U(M, fi h ,fi s ,fif) N~ 2r) s DY ^n_,/i s J . 



(65) 



Given the simple N dependence of this result, we can transform it back to momentum space 
analytically using 



1 pc+ico 

— / dNz~ N N- 2 ^ = (-ln^- 1+2 " 

2%i J c-ioQ 

and the fact that 



(1-^)1-2^ r(2r]) 



l + O 



(a-*)') 



/ M 2 \ / M 2 \ 



M 2 



(66) 



(67) 



Dropping the corrections of order (1 — z) 2 in the last step in (66), the inverted moment- 
space result can be written in a form that is identical to (50) up to an overall factor y/z, which 
amounts to a first-order power correction in the threshold region. In the limit /ih — Ms = A*/, in 
which the resummed expression reproduces the leading singular terms in the fixed-order result, 
the momentum-space formulation derived in the present work gives the singular distributions 
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Figure 10: Comparison of the moment-space hard-scattering coefficient (65) at fi^ — jj,f — M 
and fi s = M/N with M = 8 GeV (solid) with the expression (69) used in the traditional 
approach [55] (dashed). The light, medium, and dark lines correspond to LO, NLO, and 
NNLO, respectively. The dashed curves are barely visible because they are almost on top of 
the solid lines. 



in the hard-scattering kernel in precisely the form in which they appear in (9) and (73). The 
large-iV expansion in Mellin space, on the other hand, gives an expression that is obtained 
from this by the replacement 



(68) 



r m ' 




ln n (M 2 In 2 z/fij) 


l-z_ 
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— In z 



with L z = ln[M 2 (l — z) 2 /fi 2 z]. While the two expressions agree up to power-suppressed terms 
near the partonic threshold, they differ significantly for small z. This leads to large corrections 
when matching with the fixed-order results in cases where the small- z region is relevant. 
In the traditional approach to resummation, the moment-space result is written as 



C N (M 2 ,fi f ) = g (M 2 , H ) exp[G N (M 2 ,» f ))+0 



(69) 



where g$ collects the iV-independent contributions. The resummation exponent has the form 



Gn(M ,fj,f) = I dz 
'o 



1 z N-l _ x 



1-Z 



k 



— 2A (a s (k)) + D{a s (M(l-z))) 



(70) 



The coefficient A(a s ) is the cusp anomalous dimension, while D(a s ) is related to the anomalous 
dimension of the soft function. The resummation exponent evaluated to NNLO can be found 
in [55, 56, 57, 58, 59]. Note that (70) is not well defined as it stands, since the coupling constant 
is integrated over the Landau pole. The resulting ambiguity corresponds to a spurious first- 
order power correction in Aqcd/Mx> while on general grounds the cross sections is expected 
to receive power corrections starting at second order [60]. For DIS, we have shown that for 
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Figure 11: Left: Results for the Drell-Yan i^-factor obtained using threshold resummation in 
momentum-space (solid) and moment space (dashed). The default scale choices fih = 1.5M, 
fi s = (n\ + /i^)/2, and /if — M are used in both cases. Right: Moment-space results for the 
.fT-factor obtained from (65) with scale choices ^ = 1.5M and fi s = M/N (dashed) and the 
default choices \ih = 1.5M and fi s = (fj} a + /i* 1 )/^ (solid). The light, medium, and dark lines 
correspond to LO, NLO, and NNLO, respectively. 

the scale choices = M and fi s = M/N, which are implicit in the traditional scheme, the 
two approaches to resummation are equivalent up to 1/N corrections, and we have derived the 
relationship between the anomalous dimensions in the effective theory and the resummation 
coefficients in the traditional approach [22]. With the same technique, we obtain 

Did* ) p^lE^ TVl , OV7\ _ 1 

e 2 7E v r(1 + 2V) = iw{as) + v ln?DY(0; ^ _ e _ 1 + v ^ T cusp (a s ) , (71) 

where a s = a s (/x), and V = d/dliafi 2 = [P(a s )/2]d/da s . Using this relation, we reproduce 
the perturbative expression for the function D(a s ) given in [55, 56] up to third order in a^. 
In particular, for the first nonzero term we find 

D{a s ) = (2< + 2n 2 C F p ) (^f + 0(a 3 s ) . (72) 

In Figure 10, we compare our moment-space expression (65) for the hard-scattering kernel 
Cn evaluated with the scale choices fih = M and fi s = M/N with expression (69) evaluated 
as described in [55]. The results are almost indistinguishable, demonstrating that the 1/N- 
suppressed differences between the two formulations are numerically very small, and that the 
traditional approach to threshold resummation can be viewed as a special case of our RG-based 
framework. 

The fact that the moment-space resummation approach differs from the momentum-space 
approach by first-order corrections in (1 — z) leads to visible numerical effects. In order to 
illustrate this we take the inverse Mellin transformation of (63) to compute the Drell-Yan 
.fT-factor. For simplicity, we fit the polynomial form ff(y) = y a (l — y) b {l + cy + dy 2 ) to the 
parton luminosity function defined in (14), so that the product of the moments of the PDFs 
in (63) can be evaluated in closed form. The first plot in Figure 11 shows that for the same 
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choices of matching scales the moment-space approach leads to somewhat larger results for 
the .fT-factor than the momentum-space approach, and that the convergence is slightly better 
in the second case. 

In the traditional moment-space approach based on expression (70) one implicitly makes 
the scale choice fj, s = M/N, which leads to a Landau pole at N ~ M/Aq C d in the Mellin 
inversion. To see how the results obtained with this choice compare with our default results 
we adopt the so-called minimal prescription [61], which amounts to choosing an inversion 
contour in the iV-plane that does not include this pole. When this is done, we find that the 
moment-space formula (65) evaluated with ji h = 1.5M and ji s = M/N gives numerical results 
similar to those found using our default scale choices once we go beyond leading order. This 
is illustrated in the right plot in Figure 11. 

6 Conclusions 

We have presented a detailed analysis of threshold resummation for Drell-Yan production. 
Instead of the conventional Mellin moment formalism, we have performed the resummation 
directly in momentum space, using RG evolution. We have obtained analytic expressions for 
the resummed hard-scattering kernels that are free of the Landau-pole ambiguities inherent 
in the traditional approach. After deriving the necessary matching coefficients and anomalous 
dimensions from known perturbative results in the literature, we have performed the resum- 
mation for the cross section and rapidity distribution at NNLO in RG-improved perturbation 
theory, corresponding to N 3 LL accuracy. 

The main goal of this work was to address the question to what extent resummation can be 
phenomenologically relevant given that in all practical applications the true threshold region 
is experimentally not accessible, because the parton distributions fall off very steeply at large 
x. It has been argued in the literature that precisely this fall-off forces the Drell-Yan process 
to the partonic threshold region, in which large logarithmic corrections arise. Our analysis 
confirms the existence of this effect and quantifies its importance. In the true endpoint region, 
we find that the scale of soft radiation is an order of magnitude smaller than the naive estimate. 
The extra suppression factor is given by the sum of the exponents of the fall-off of the quark 
and anti-quark distributions near x — 1. On the other hand, for very small lepton-pair masses 
this effect becomes inoperative because the fall-off of the parton distributions is not very steep 
at low x. An analysis of the convergence of perturbation theory in the intermediate region 
shows that resummation effects become relevant for M larger than about 0.4-^/s. At the 
largest measured Drell-Yan masses resummation effects can thus be significant. In the case of 
the experiment E866/NuSea, a fixed-target experiment with pp collisions at y/s = 38.76 GeV, 
these effects increase the NNLO (NLO) cross section by 7% (25%) at M = 16 GeV. 

We do not confirm the recent claim of large negative resummation effects for the rapidity 
distribution measured in the E866/NuSea experiment at M — 8 GeV [15]. We find that the 
fixed-order NNLO threshold contribution is positive and a good approximation to the full 
NNLO correction. Higher-order terms beyond NNLO turn out to be negligible, so that it 
is unnecessary to resum them. To check whether the large effect seen in [15] could be an 
artifact of the moment-space formalism, we have applied our approach also in moment space 
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and showed that it encompasses the traditional moment-space approach. We find that the two 
methods give compatible results, and that in both schemes resummation increases the cross 
section. The Drell-Yan cross section is dominated by the contribution of the threshold region 
even at small masses of the lepton pair. For example, for M = 8 GeV the threshold terms give 
rise to 93% of the NLO correction to the total cross section and 97% of the NNLO contribution. 
If this is a more generic feature of hard cross sections, then threshold resummation techniques 
provide an efficient way to obtain the dominant part of the higher-order corrections. 

Perhaps the most important outcome of our analysis is a quantitative understanding of 
the emergence of an effective physical scale characterizing the soft radiation in the Drell-Yan 
process. This scale is generated through an intricate interplay of dynamical and kinematical 
effects. It will be interesting to explore how the same mechanism affects other collider processes 
such as Higgs production at the Tevatron and LHC Furthermore, it would be useful to analyze 
resummation for tt-production, which is phenomenologically relevant, since the Tevatron now 
produces top-quark pairs with quite high invariant masses. There was some controversy on 
how to best perform the resummmation in this case [70, 71, 61, 72]. It would be interesting 
to revisit the issue using effective field-theory methods. 
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A Two-loop coefficients of the leading singular terms 

The leading singular terms in the fixed-order perturbative expressions for the hard-scattering 
kernels C^- can be derived by using the perturbative expansions of the hard and soft functions 
in (15). In this way, we obtain for the coefficients Pi{z) in (9) 
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9 r9 / 93 3tt 2 
~ L + <-16 + - 



. 511 33tt 2 23vr 4 15Cs 

+ 5 ^) L + -^-^ + 120"- — 



1-z 

11 



+ — L + 3L — 8 — 



2tT 



1-2 



+ 



--L 2 + 
16 



193 llvr 2 



48 



36 



1-z 
.2 



3C 3 \ L 1535 47vr 2 23tt 4 43Cs 
~2~ ) ~792~ + l36 720" + T2 _ 



+ 



101 llvr 2 7C 3 

1 h — 

27 18 2 



1 


/67 tt 2 \ 

/U-5-J 


" £ 2 " 


11 


[ L ' " 




.1-2. 


.1-2. 




.1-2. 


+ 



P^) =5(1-2) 



1 2 ( 17 vr 2 \ r 127 4vr 2 Cs 

-L 2 + 1 LH h — 

4 V 12 9 } 48 9 3 



'28 _ 27T 2 
1 27 ~ T" 



1 - z 



10 



1-2 



1 

+ 3 



1-2 



(73) 



As before, we use the abbreviations L = ln(M 2 //x 2 ) and L z = ln[M 2 (l — z) 2 /fi"jz]. If desired, 
the plus distributions involving powers of L z can be reduced to distributions of the form 
[ln n (l - z)/(\ - z)] + using the identity 

[/(2) g(z)} + = f(z) [g(z)} + - 5(1 - z) f dz> f(z') [g(z')} + . (74) 

Jo 

B Matching coefficients and anomalous dimensions 



For completeness we list the perturbative expansions of the various matching coefficients and 
anomalous dimensions required to evaluate our RG-improved result (50) at NNLO. 



B.l Two- loop matching coefficients 

The matching conditions for the Wilson coefficient Cy evaluated at time-like momentum 
transfer can be obtained by analytic continuation from the corresponding expression valid at 
space-like momentum transfer. Using the known two-loop result for the on-shell QCD form 
factor [62, 63, 64, 65], we find [22] 

CV(-M 2 -*e,/i) = 1 + ^ (-L 2 + 3L - 8 + ^ +C F {^f [C F H F + C A H A + T F n f H f ] , 

(75) 
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where 

L 4 ors /^25 tt 2 \ r2 / 45 3tt 2 o ^ \ r 255 7tt 2 83tt 4 o ^ 

H F = 3L 3 + L 2 + + 24C 3 L + + 30C 3 , 

F 2 V 2 6 / V 2 2 s 7 8 2 360 S3 ' 

11 r3 ( 233 ^ r2 ^545 llvT 2 _ \ r 

51157 337?r 2 IItt 4 313 
~ ~648 108~ + + IT ^ 3 ' 
4 „ 38 r2 / 418 4tt 2 \ r 4085 23tt 2 4 > 

and L = ln(M 2 /yU 2 ) — ?7r. This result agrees with the corresponding expression given in [36]. 

The matching condition for the function Sdy can be derived most easily using relation 
(48), which relates it to the perturbative expansion of the position-space Wilson loop Wdy at 
time-like separation. We can then use the explicit two-loop expression for the position-space 
Wilson loop obtained in [39] . To this end, we have obtained the fourth-order coefficient in the 
expansion of a certain Appell hypergeometric function, 



The resulting two-loop expression for the soft function reads 



s DY (L, /i) = 1 + ^ (2L 2 + + Cf (^) 2 [C F W F + C A W A + T F n f W f ] , 



(77) 



(78) 



where 



F 3 18 2\ 3 

22 ^ 3 /134 2tt 2 \ r0 ( 808 . _ \ . 2428 67tt 2 tt 4 22 



W A = L 6 + £ + + 28C3 L + + Cs , 

9 V 9 3 J V 27 / 81 54 3 9 S3 ' 

8 r3 40 r2 224 r 656 IOtt 2 8 . 

W f = - Lr L H L h - G . (79) 

/ 9 9 27 81 27 9 s y ' 

Note that the C F terms exponentiate, which is a consequence of the non-abelian exponentiation 
theorem for Wilson loops [66, 67]. Our result for the function s DY agrees with a corresponding 
expression entering in the moment-space resummation approach studied in 



B.2 Three-loop anomalous dimensions 

Here we list expressions for the anomalous dimensions and the QCD /3-function, quoting 
all results in the MS renormalization scheme. We define the expansion coefficients of the 
anomalous dimensions and the QCD /3-function as 

^ / % „ a s „ { ot s \ 2 _ fa s \ 3 
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and similarly for the other anomalous dimensions. 

The expansion of the cusp anomalous dimension r cusp to two-loop order was obtained 
some time ago [46], while recently the three- loop coefficient has been calculated in [49]. For 
the four-loop coefficient T 3 we use the Pade approximant T 3 = r 2 /Ti. The results are 



r = 4cv, 
r\ = ac f 

To = 4CV 



67 tt 2 \ 20 



cl 



245 



134tt 2 IItt 4 22 . ., ,. 



418 40tt 2 _ 56 

"27" + ~W ~ IP 3 



55 



+ C F T F n f [-— + l6( 3 



27 T ^ 



T 3 « 7849, 4313, 1553 for n f = 3, 4, 5 . 



(81) 



The anomalous dimension 7 y can be determined up to three-loop order from the partial three- 
loop expression for the on-shell quark form factor in QCD, which has recently been obtained 
in [65]. We find 
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(82) 



V 729 27 27 

The anomalous dimension 7^ is know to three-loop order from the NNLO calculation of the 
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Altarelli-Parisi splitting functions [49]. The expansion coefficients are 
lt = 3C F: 

it = Cl (1 - 2^ + 24(a) + C F C A (j + 2 ^f- 12 ^) - C ^ n f (I + IT 

7^^(| + 3^ + ^ + 68 C 3-^C3 - 240C 5 

/151 205jt 2 247?r 4 844 , 8?r 2 , ,„„, 

+ c ^{- - — - -nr + - 3-6 + -6 + 120& 

F ; V 9 135 3 s 

n ^2 f !657 2248tt 2 tt 4 1552 , 

^ ^ _ / 1336tt 2 2vr 4 400 A 
l , t2 2 ( 68 1607T 2 64 \ 

Using these results, one can compute the expansion coefficients for the anomalous dimension 
7 W of the Drell-Yan soft function from the relation ^ w = 1^ + 7 y . This yields 

/224 4tt 2 



lY = C F C A (-^ + ^ + 28C 3 ) + C F T F n f ( 



27 



(84) 



We do not list the three-loop coefficient. 

Finally, the expansion coefficients for the QCD /5-function to four-loop order are 

11 4 
A) = y C A - - T F n f , 

34 20 

A = - - - C A T F rif - 4CV7W , (85) 



3 " 3 

h Al + (2Cl-^C F C> 1415 



54 A V * 9 27 



C^) T F rif +(y Cf+1 YT Ca ) T ^ ' 

149753 , / 1078361 6508 , \ /50065 6472 > \ 2 1093 

A = — + 3564 Cs - + _ Cs j n/ + (_ + _ C3 j „} + _ : 

The value of /3 3 is taken from [68] and corresponds to N c = 3 and Tp = \. 
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B.3 Renormalizat ion-group functions 

We now give the perturbative expansions of the functions S and ar defined in (40), working 
consistently at NNLO in RG-improved perturbation theory. At this order we need to keep 
terms through 0{at 2 s ) in the final expressions. The resulting expression for ar is given by 



r I a s (fx) fT 1 ft\ a s (fx) - a s (u) 



+ 



Ei _ ft _ h (Ei 
r ft ft Vr 



ft 

ft 



32tt 2 



+ 



(86) 



Similar expressions with the Tj replaced by the coefficients 7^ or 7* hold for the functions a 7 v 
and a 7 ^, respectively. The NNLO expression for the Sudakov exponent S is more complicated. 
It reads [22] 
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(87) 



where r = a s {Li)/a s {v). Whereas the three- loop anomalous dimensions and ftfunction are 
required in (86), the expression for S also involves the four-loop coefficients T 3 and ft. 
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Figure 12: Example of a cut diagram contribution to the Drell-Yan process. 

C Cut diagrams in the Keldysh formalism 

To perform the effective theory analysis of the process, we would like to have a path-integral 
definition of the hadronic quantity of interest. In the case of DIS (and also for inclusive B- 
decays), one considers the discontinuity of forward matrix elements of time-ordered products 
of the electroweak currents and studies their factorization properties. The Drell-Yan cross 
section cannot be written in a similar form. The reason is that not all cuts of the relevant 
Feynman graphs correspond to the same physical process. For example, in addition to the 
contribution from the cut indicated in Figure 12, the discontinuity of the same diagram also 
gets a contribution from a cut through the triangle loop on the left, which describes p+p — > X 
with a virtual lepton pair. 

A path-integral method for the direct evaluation of cut diagrams is the Keldysh (or time- 
loop) formalism [37, 38]. Let us illustrate the method for a scalar field theory with field 4>(x). 
Instead of the usual action, where one integrates over time from t = —oo . . . oo, one considers 
a path integral with an action 

5(0)= JdtCtf), (88) 

where the contour C first runs from t = — oo + i5 ... oo + i5 and then back, with a negative 
imaginary part, along the contour shown in Figure 13. When evaluating expectation values 
of fields, one needs to specify on which part of the contour the fields reside. We denote the 
fields living on the first half of the contour by (f)+(x) = <p(x + i5) and the ones on the second 
half by (f>-(x) = (f>(x — i5). The path- integral expectation values of fields correspond to path- 
ordered vacuum expectation values. Since the fields <f> + are positioned at earlier points along 
the contour, they are to the right of all fields 0_. Furthermore path-ordering translates into 
anti-time-ordering on the second part of the contour, so that we find 

J V(f>4> + (xi) . . . 4>+{x n ) (f)-{xi) . . . <f>-(x m ) exp {iS(4>)} 

= (0\T{<f>(x 1 )...<f>(x m )}T{ ( f ) (x 1 )... ( f>(x n )}\0). (89) 
Here T denotes anti-time-ordering. Note that under a field redefinition <j)(x) — > f((f>(x)) (f>(x), 
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Figure 13: Contour in the complex t-plane along which the Lagrangian is integrated to ob- 
tain the action in the Keldysh formalism. (To dampen the oscillatory behavior of the path 
integrand, the entire contour should further be rotated clockwise by a small angle.) 

the plus and minus fields transform separately <f>±(x) — > f{<j)±{x)) <f>±(x). This is relevant 
when one uses a field redefinition to decouple the soft gluon fields from the quark fields in the 
current operator. 

The Keldysh formalism is useful, because it gives a path integral formulation of squared 
amplitudes, e.g. 

J2(^(Px - P) I (X\ T . . . cf>(x n )} |0) | 2 

x 

= J d A x e~ ip - x (0| T {0(xi +x)... <j)(x n + x)} T {<j>( Xl ) . . . <j)(x n )} |0> , (90) 

which is then rewritten as an expectation value of <f>± fields using (89). 

For our perturbative analysis, we need the Feynman rules to calculate matrix elements of 
the form (89). The rules are simple: those for the field <p + are the usual Feynman rules and 
those for 0_ are the complex conjugate of the usual rules. Since all interactions are local, there 
are no vertices involving both 0+ and 0_ fields. The only connection is the cut propagator 

<0|<M*)<M0)|0> = j -0^e-^(27r)5(p 2 -m 2 )6(p°). (91) 

This expression is familiar from the Cutkosky rules [69] used to extract the contribution of a 
given cut to an ordinary Feynman diagram. 

Considering the case of the Drell-Yan cross section in (18), the relevant quantity to consider 
in the Keldysh formalism is the current product J^(x) J +fl (0), where = ^ q e q q±^q±. 
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